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SECTION  1 


INTRODUCTION 

It  is  generally  recognized,  on  both  theoretical  grounds  and 
evidence  from  the  bombing  of  Hiroshima  (Bond,  1946),  that  nuclear 
weapons  might  induce  simultaneous  ignitions  over  very  large  areas,  which 
could  develop  into  very  large  fires.  Such  fires  create  huge  central 
convective  plumes  which  rise  to  large  heights.  It  is  also  believed  that 
such  fires  create  very  strong  radially-inward  winds  in  the  atmospheric 
surface  layer  which  increase  mixing  and  provide  fresh  air  to  the  burning 
area  but  inhibit  radially  outward  spread  of  the  fire  (Carrier,  et  al., 
1985). 

Such  large  fires  could  have  a  profound  effect  on  the  earth’s 
atmosphere.  Smoke  particles  injected  into  the  stratosphere  could  widely 
spread,  blocking  solar  radiation  and  thus  causing  a  considerable  drop  in 
temperature  at  the  earth's  surface.  Attempts  to  quantitatively  predict 
the  effect  of  such  firestorms  are  greatly  impeded  by  many  uncertainties 
(Carrier,  1985),  some  of  which  are  related  to  the  characteristics  of 
large-area  fire  plumes. 

The  purpose  of  the  present  study  is  to  investigate  features  of  the 
fire-induced  flow  by  physical  simulations  in  small-scale  models. 
Current  thoughts  on  the  gross  characteristics  of  large-area  fires,  based 
primarily  on  previous  studies  of  related  phenomena,  are  briefly 
described  in  the  beginning  of  the  report.  The  rationale  of  the  physical 
simulation  schemes  used  in  this  investigation  and  the  dimensionless 
description  of  the  flow  are  then  discussed.  Finally,  the  experimental 
results,  which  depict  important  properties  of  large-area  fires  in  a  calm 
atmosphere  are  presented  and  discussed. 


SECTION  2 


GROSS  FEATURES  OF  THE  FIRE- INDUCED  FLOW  FIELD 


A  large-area  fire  (LAF)  is  defined  in  this  work  as  one  which  burns 
simultaneously  over  a  large  area  of  many  square  kilometers  and  in  which 
the  height  of  the  average  combustion  layer  is  small  compared  to  its 
horizontal  size.  The  correlation  of  Thomas  (1963),  shown  in  Fig.  2,  can 
be  used  to  determine  whether  a  given  fire  can  be  classified  as  a  LAF 
according  to  this  definition. 

Figure  1  schematically  describes  our  assumed  structure  of  an 
axisymmetric  LAF  in  a  calm  atmosphere.  In  the  figure  we  have  identified 
four  zones  of  distinct  characteristics.  The  shallow  combustion  layer 
(CL)  is  the  immediate  layer  near  the  ground  which  is  limited  to  the 
average  flame  height.  The  region  above  the  combustion  layer  is  divided 
into  two  zones:  the  central  core ,  which  covers  approximately  one-fifth 
of  the  burning  area,  and  the  radial  ( convective)  boundary  layer ,  where 
radially-inward  winds  carry  heat,  smoke,  as  well  as  fresh  air,  toward 


the  central  core.  Above  the  central  core,  at  z  >  Z  ,  where  Z  is  the 

P  P 


pinch  point,  a  convective  plume  is  formed,  which  could  rise  to  large 
heights . 

It  is  noted  that,  in  contrast  to  the  above  schematic  shape  of  a 
large-area  fire  plume,  Small  et  al.  (1980,  1984)  showed  a  relatively 
small  necking  of  the  rising  plume.  Carrier  et  al.  (1984)  showed 
theoretically  that  considerable  necking  of  the  plume  must  occur, 
however,  the  relative  height  of  the  combustion  layer  in  some  schematic 
descriptions  of  large  fire  plumes  shown  in  their  work  is  large. 

The  nature  of  the  convective  plume  at  large  heights  above  the 
ground  is  expected  to  be  similar  to  that  of  weak  convective  plumes.  The 
behavior  of  such  plumes,  in  either  neutrally  stratified  or  in  stably 
stratified  atmospheres,  has  already  been  treated  by  many  investigators 
(Rouse,  et  al.,  1952;  Priestley  and  Ball,  1955;  Morton,  1965;  Morton, 
et  al.,  1956;  Briggs,  1969;  and  others).  These  studies  are  conveniently 
summarized  in  Turner  (1973). 

The  major  parameter  governing  the  flow  of  a  weak  convective  plume 


is  the  total  buoyancy  flux  Fq,  which  remains  constant  in  a  laboratory 


environment  of  constant  density  or  in  an  ideal,  adiabatic  atmosphere. 
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Vs, 


»vv 


Assuming  a  self-similar  solution,  one  obtains,  for  this  case,  (Turner, 
1973)  that 


w  =  f  (r) 

~l/3  —1/3  rl  V  ’ 

r  z 

(i) 

0 

_ _  =  f 

P2/3~ 5/3  r2  V 

i?  z 

(2) 

and 

b/z  =  P  (3) 

where  z  is  the  height  from  the  origin  of  the  plume  and  b  is  a 
characteristic  size  (radius)  of  the  plume.  The  linear  growth  of  b 
implies  that  the  mean  inflow  velocity  across  the  edge  of  the  plume  is 
proportional  to  the  characteristic  local  mean  velocity  (w)  at  each 
section  (Batchelor,  1954)  or: 

(b2w)  =  2abw  .  (4) 

Values  of  a  (=  6p/5)  ranging  from  0.093  to  0.12  have  been  suggested  in 
the  literature. 

Experiments  indicate  that  both  f^  and  f^  are  approximately 
Gaussian.  Using  the  data  of  Rouse  et  al.  (1952),  for  a  point  source  at 
ground  level,  they  can  be  written  as: 

2 

f.  =  4.7  exp  ( - - - y)  (5) 

2(0.072z r 


f9  =  11  exp(-  - - - x)  =  11  exp(-  In  2  — )  (6) 

2(0 . 084z;  bj  5Q 

where  b^  =  O.lz  is  the  radius  where  g'(r)/g'(0)  =  0.5. 

The  above  formulation  assumes  an  ideal  point  source  of  buoyancy 
flux  near  the  ground;  namely  a  small  initial  value  of  Ap/p  ,  a  small 
initial  mass  flux  and  zero  initial  momentum  flux.  However,  it  is  also 
expected  to  be  valid,  at  large  heights,  for  finite  size  sources  with 
initially  large  values  of  Ap/p  ,  mass  flux  and  momentum  flux,  except 

3 


that  the  self-similar  plume  at  large  heights  would  appear  to  have 
emerged  from  a  virtual  origin  at  z  =  Zv,  which  is  not  identical  to  the 
physical  origin  of  the  source,  as  shown  in  Fig.  1.  The  value  of  z  in 
Eqs.  ( 1 ) - (6)  is  thus  z  =  z  -  Zy.  Large  initial  plume  size  and  an 
initial  momentum  flux  lower  the  height  of  the  virtual  origin  of  the 
plume,  whereas  large  values  of  Ap/p  at  the  source,  which  attenuate 
the  initial  mixing,  raise  the  height  of  the  virtual  origin.  It  is  thus 
expected  that  Zy  for  large* area  sources  would  be  negative. 

An  analytical  solution  describing  a  plume  from  a  circular  area 
source  of  heat  was  proposed  by  Priestley  and  Ball  (1955).  In  the  far 
field,  the  solution  is  identical  to  that  of  Morton  et  al.  (1956),  except 
that  the  virtual  origin  of  the  far  field  plume  is  at  Zy  =  -10R.  This 
unrealistic  result  is  due  to  the  assumption  that  the  similarity  of 
solution  extends  to  z  =  0,  which  appriori  eliminates  a  possibility  for 
plume  pinching. 

Plumes  in  stable  atmospheres  lose  their  buoyancy  flux  as  the  plume 
rises  to  regions  with  smaller  potential  temperatures.  The  behavior  of 
plumes  in  a  constant  stable  density  gradient  in  the  atmosphere,  have 
been  treated  by  Morton  et  al.  (1956),  who  described  the  dimensionless 
width  (b'),  velocity  ,  acceleration  ,  and  buoyancy  flux, 

f(Zj)  =  Fo(z1)/Fo(o),  (7) 


as  a  function  of  a  dimensionless  height, 


Z1  = 


1 


0.41  a  "1/2  F  1/4  N"3/4 
G  o 


(8) 


As  pointed  out  by  Morton  et  al.,  the  initial  plume-spread,  up  to 
approximately  z^  =  2,  is  linear  and  the  solutions  are  little  different 
from  those  in  neutral  environments.  At  approximately  z^  =  2.8,  however, 
the  vertical  velocity  vanishes  and  a  radial  spread  of  the  plume  starts. 
It  is  also  interesting  to  note  that  the  initial  loss  of  buoyancy  flux  of 
the  plume  is  very  small  and  most  of  the  buoyancy  flux  is  lost  after  the 
width  of  the  plume  becomes  large.  For  example,  f(0.5)  =  0.98,  f(l)  = 
0.86,  f (1 .5)  =  0.60,  and  f (2 .0)  =  0.04. 


4 


Briggs  (1969)  collected  "plume-rise"  data  from  point  sources  in  a 
calm  stable  atmosphere  and  showed  that 


f1/4  g-3/8 


3  5.0. 


The  value  of  F  =  FQ(D)/7t  is  used  in  the  above  equation,  and  many  other 


plume-rise  equations,  so  that  Briggs  plume  rise  corresponds  to  approxi¬ 


mately  Zj  =  2.5. 


It  has  also  been  established  that  atmospheric  circulation  could 
have  a  profound  effect  on  convective  plumes  (Long,  1966;  Emmons  and 
Ying,  1967,  Carrier  et  al.,  1985).  The  radially- inward  flux  of  wind 
will  enhance  the  circular  atmospheric  motion  to  preserve  angular  momen¬ 
tum.  This  circulatory  motion  could  reduce  the  entrainment  of  ambient 
air  into  the  lighter  core  of  the  convective  plume,  reduce  its  spread  and 
considerably  increase  its  rise  in  a  stably  stratified  atmosphere. 

Radiative  cooling  of  the  hot  convective  plume,  which  depends  on  its 
temperature  and  smoke  contents,  can  also  reduce  its  buoyant  flux  and 
decreases  the  plume  rise.  Radiative  cooling  is  not  included  in  this 
study  and  the  reader  is  referred  to  the  extensive  review  by  Kanury 
(1984).  Based  on  this  review  it  is  estimated  that  radiative  cooling  will 
not  have  a  considerable  effect  on  the  structure  of  the  plumes  from  LAF. 
Since  the  magnitude  of  radiative  cooling  depends  on  the  temperatures  in 
the  rising  plume,  it  is  hoped  that  the  result  of  this  study  will  also  be 
useful  for  better  estimates  of  the  effect  as  well  as  for  estimating  the 
soot  fraction  in  the  smoke  produced  by  fires,  which  depends  critically 
on  whether  the  fires  are  well  ventilated  (Rathjens  and  Siegel,  1985). 
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SECTION  3 


PHYSICAL  MODELING  OF  THE  FIRE- INDUCED  FLOW  FIELD 

The  problem  of  simulating  large-scale  fire  plumes  in  small-scale 
physical  models  has  received  previous  attention  (Nielsen,  1965;  Long, 
1966;  Parker  et  al.  (1968);  Williams,  1969  and  others).  A  recent  com¬ 
prehensive  review  of  the  problems  was  prepared  by  Kanury  (1984). 
Clearly,  it  is  impossible  to  correctly  simulate  in  small  models  all  the 
physical  phenomena  which  take  place  in  large  fires  and  affect  the  com¬ 
bustion  and  flow  processes.  One  must  therefore  resort  to  approximate 
models  in  which  only  the  most  important  processes  are  simulated.  Pre¬ 
vious  experience  has  shown  that  such  approximate  physical  models  are 
successful  in  simulating  complex  three-dimensional  atmospheric  flows 
with  complex  boundary  conditions,  which  cannot  easily  be  studied  analy¬ 
tically  or  numerically  (Cermak,  1971,  1981).  Moreover,  physical  models 
provide  a  data  base  which  can  be  used  to  verify  and  improve  semi- 
empirical  numerical  and  analytical  models.  In  fact,  most  coefficients 
in  turbulent  closure  models,  as  well  as  coefficients  in  the  classical 
analytical  treatments  of  weak  plumes  and  jets,  have  been  determined  by 
small-scale  simulations. 

'  Nielsen  (1965)  and  Long  (1966)  proposed  to  produce  in  models  of 
large  fires  a  constant  ground  temperature,  T  +  AT,  where  T  is  the 
ambient  temperature.  Thus,  the  scaling  of  the  various  parameters  in  the 
model  and  the  prototype  become  a  function  of  the  ratio  of  AT  in  the 
model  to  that  outside  the  combustion  layer  in  the  prototype. 

The  simulation  proposed  by  the  principal  investigators  is  based  on 
the  observation  that  the  primary  driving  force  that  shapes  the  structure 
of  the  flow,  above  the  thin  combustion  layer,  is  the  buoyancy  flux  per 
unit  area  from  that  layer,  which  will  be  denoted  by  ct.  Assuming  that 
a  is  a  constant,  it  is  generated  in  the  models  by  a  flow  of  lighter- 
than-air  gas,  such  as  helium,  through  a  porous  plate  at  z  =  0.  It  is 
noted  that  helium  has  a  relative  density  deficiency  Ap/p  =  0.85,  which 
is  close  to  that  of  burned  fuel  at  stoichiometric  proportions  before  it 
mixes  with  additional  air  (Carrier  et  al.  1984).  Accordingly,  the 
scaling  of  all  variables  in  the  model  is  determined  by  the  geometric 


scaling  of  the  model 


\  =  VLp 

and  the  scaling  of  the  buoyancy  flux 

\  =  a  /a 
a  m  p 


(10) 


(11) 


where  \  =  x  /x  is  the  ratio  of  any  variable  x  in  the  model  (m)  and 
x  m  p 

in  the  prototype  (p). 

The  scalings  of  velocities  (U),  effective  gravitational 
accelerations  (g'),  relative  density  deficiencies  (Ap/p  )  are  therefore 


X  _  \ 1/3  ,1/3 


2/3  ,-1/3 


\  \  =  \  =  \  J  A_ 

g'  Ap/pa  Ap  a  1, 


and 


V  ■  V\>  =  t'3  Xal/3 


(12) 

(13) 


(14) 


The  present  experiments  were  restricted  to  models  of  large-area  fire 
plumes  in  adiabatic  and  calm  atmospheres,  which  are  represented  in  a 
small-scale  model  by  a  constant  density  environment.  Simulation  of 
large-area  fires  in  stably  stratified  atmospheres,  is  discussed  in  Poreh 
et  al.  (1986) . 

A  very  important  additional  requirement  for  correct  simulation  is 
that  the  flow  in  the  model  be  turbulent,  as  in  large-scale  fires.  This 
requirement  implies  that  the  Reynolds  number  of  the  model,  say 


Re  = 


a1/3  R4'3 


(15) 


be  sufficiently  large.  It  was  experimentally  found  that  when  Re,  as 
defined  above,  is  larger  than  6000,  the  flow  over  the  buoyancy  source  in 
the  model  appears  to  be  fully  turbulent. 

It  should  be  stressed  that  the  scaling  criteria  for  the  proposed 
simulation  do  not  in  any  way  contradict  those  proposed  by  Long,  except 


that  the  condition 


is  not  considered  by  the  authors  to  be  an  absolute  necessary  condition 
for  modelling.  Such  a  relaxation  has  one  important  advantage.  It 
enables  one  to  increase  the  buoyancy  flux  in  the  model  and  thus  increase 
both  the  induced  velocities  and  the  Reynolds  number  of  the  model. 
However,  it  is  not  acceptable  for  cases  where  strong  swirling  of  the 
plume  occurs.  To  model  correctly  the  effect  of  ambient  air  vorticity, 
its  value  in  the  model  should  scale  as 

\i  ■  \//3  V2/3  •  (17) 

In  addition,  it  is  essential  that  the  centripetal  accelerations  acting 

2 

on  lighter-than-air  fluid  ApV./(p  r),  where  V.  is  the  circumferential 

*P  a  <p  2 

velocity,  be  scaled  in  the  model  as  g'  (Eq.  (13)).  Since  VT7r  scales 

2/3  -1/3  * 

as  a  R  ,  such  models  must  satisfy  Eq.  (16)  or 

\  *  *a  •  <18) 


This  requirement  ensures  that  the  radial  densimetric  Froude  number, 

associated  with  the  centripetal  acceleration,  is  the  same  in  the  model 

and  prototype.  The  densimetric  Froude  number  for  flows  with  horizontal 

1/2 

velocities  is  usually  defined  as  U/(g'L)  .  In  our  case  the  mean 

2 

vertical  velocity  w  replaces  U,  ApV./(R  p  )  replaces  g'  and  R  replaces 

™  a  2  1/2 
L,  so  that  the  densimetric  Froude  number  becomes  w/fV^  Ap/pa) 

This  dimensionless  number  becomes  important  only  when  its  value  is 

small.  When  a  strong  swirling  motion  is  not  present  and  <  w,  the 

effect  of  Ap/p  ,  which  is  usually  <  0.3,  is  small,  and  there  is  no  need 

to  satisfy  Eqs.  (17)  and  (18)  in  the  model. 

One  may  also  generate  buoyancy  flux  from  area  sources  using  hot 

plates  with  a  constant  power  supply  Q^.  Approximately,  the  buoyancy 

flux  per  unit  area  in  such  models  would  be  Qrg/(p  c  T  A),,  where  A  is  the 

L  3  p  3 

area  of  the  source.  Such  a  1:537  scale-model  was  used  by  Parker  et  al. 
(1968)  to  simulate  the  335  m  x  320  m  Flambeau  Fire.  The  primary  dis¬ 
advantage  of  such  models  is  that  an  unknown  part  of  the  heat  flux  is 
lost  by  radiation.  Measurements  by  Lee  (1969)  showed  that  radiation 
losses  can  be  as  large  as  42  percent. 


The  use  of  small  fires  to  simulate  large  area  fires  is  also 
subjected  to  the  above  limitation.  In  addition,  the  relative  thickness 
of  the  combustion  layer  in  small-scale  fires  is  usually  larger  than  in 
large-area  fires  and  thus  some  of  the  buoyancy  in  such  models  is  pro¬ 
duced  at  relatively  large  heights.  This  would  probably  raise  the  loca¬ 
tion  of  the  virtual  origin  of  the  far  field  plume  and  effect  the  magni¬ 
tude  of  the  induced  horizontal  velocities  near  the  ground.  One  may 
decrease  the  thickness  of  the  combustion  layer  in  such  models  by  reduc¬ 
ing  the  burning  rate,  however,  the  induced  horizontal  velocities  might, 
in  such  cases,  be  very  small. 

In  the  proposed  model,  the  entire  burning  process  is  absent,  and 
buoyancy  is  produced  at  ground  level  by  injection  of  helium,  which 
roughly  simulates  burned  fuel  before  it  mixes  with  additional  air.  One 
may,  however,  define  a  pseudo-combustion  layer  (PCL)  in  such  models  by 
the  distance  from  the  ground  where  the  helium  is  mixed  with  about  six 
times  air  by  mass,  as  this  is  usually  a  rough  criterion  for  the  thick¬ 
ness  of  the  combustion  layer  in  fires  (Thomas,  1963;  Carrier  et  al., 
1984). 

By  measuring' the  thickness  of  the  PCL  in  the  model,  the  validity  of 
such  simulations  can  be  examined.  Since  it  is  assumed  that  the  relative 
height  of  the  combustion  layer  in  LAF  is  small,  it  is  absolutely  neces¬ 
sary  that  the  PCL  in  the  model  be  relatively  small.  This  requirement 
could  set  an  upper  limit  to  the  value  of  o  in  the  model.  As  will  be 
shown  later,  the  models  in  the  present  investigation  satisfy  this 
requirement. 

Our  basic  assumption,  that  the  flow  field  above  the  combustion 
layer  is  determined  by  the  buoyancy  flux  per  unit  area  and  the  size  of 
the  fire,  implying  that  all  the  velocities  in  the  fire-induced  flow 
field,  in  both  the  model  and  the  prototype,  are  proportional  to 
(ctR)1/3,  or 


u  (r,z)  _  ,,r  z. 
(aR)1'3  "  ' 


(19) 


The  equivalent  dimensionless  equation,  which  uses  the  total  buoyancy 


flux  of  the  fire  F  =  nR  o  and  R  rather  than  o  and  R,  is 
o 


u  (r,  z) 


=  t(§  ;  I)  • 


n-l/3  F 1 / 3  r-1/3  lvR  ’  R 
o 


(20) 


The  flow  in  the  convective  plume  is  expected  to  be  a  function  of  the 
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total  buoyancy  flux  Fq  =  TiR  a  and  independent  of  a  and  R.  Thus,  at 
this  region 


f(i ;  =  (i)_1/3  F(^ 

z 

which  is  consistent  with  Eq.  (1).  Similarly, 


_ S _  =  _ S _  = 

(no)2'3  R-1/3  F2/3  R'5/3  8(r  ’  R) 

o 


In  the  convective  plume  Eq.  (22)  reduces  to  Eq.  (2). 

In  models  where  Mg')  =  \(Ap/p  )  =  1  and  Eq.  (17)  is  satisfied,  the 

3  1/2 

induced  velocities  in  the  model  would  be  proportional  to  R  as 
suggested  by  Long  (1966)  and  Parker  et  al.  (1968). 

The  flow  above  the  combustion  layer  in  the  outer  region  of  the 
radial  convective  boundary  layer,  where  r  =  0(R),  is  expected  to  be 

determined  by  the  distance  x  from  the  outer  edge  of  the  fire,  by  a 

1/3 

and  by* the  approach  velocity  UR,  which  is  proportional  to  (aR) 
Thus,  in  this  region, 

L  r  f  (?  ££_)  (23) 

UD  Vx  »  17  3J  ’ 

K  UR 


^  =  Ml  >  •  (24) 

R 

If  one  assumes  that  the  flow  in  this  region  is  self-similar,  in  the 


sense  that 


—  =  ?(-?—) 
UR  US(x)J 


UR  _  z 
a  '  8  6(x)  ’ 


the  characteristic  thickness  of  the  convective  layer  would  be 


I*' 


r*,-, 
ft 1 


Br 

$ 


' 

& 

:«• 
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HIM! 


Since  is  relatively  small,  one  expects  6  to  be  proportional  to 

Ur*1,  or 

=  .  (28) 

R 

When  both  ambient  winds  (U  )  and  vorticity  (ft  )  exist  in  the  area, 

a  a 

they  should  be  scaled  in  the  model  according  to  Eqs.  (12)  and  (17)  and 
the  dimensionless  description  of  the  flow  would  be: 

/  u 
u  _  f  /  r  z  a 

ffl/3  R1'3  ’  »  ’  a1/3  R1'3 

It  is  estimated  that  the  dimensionless  ambient  vorticies  in  the 

laboratory  experiments  of  Emmons  and  Ying  (1967)  were  one  order  of 

magnitude  larger  than  those  usually  found  in  the  atmosphere  over  regions 

which  are  larger  than  the  considered  LAF.  Thus,,  the  drastic  effects  of 

rotation  found  in  that  study  are  generally  not  expected  to  exist  in  most 

LAF  plumes  in  neutrally  stable  atmospheres. 

It  is  also  interesting  to  note  that  the  flow  field  generated  by  a 

heat  island  is  similar,  within  the  approximation  of  our  model,  to  the 

flow  field  generated  by  large-area  fires  with  very  shallow  combustion 

1/3  1/3 

layers,  except  that  U  /a  R  for  a  heat  island  is  usually  many 

times  larger  than  in  large  fires,  and  usually  cannot  be  ignored.  Physi¬ 
cal  simulations  of  heat  islands  (Yamada  and  Meroney,  1971;  Faust,  1981) 
are  therefore  approximate  models  of  large  fires  and  their  results  can  be 
used  in  the  study  of  large  fires,  provided  that  the  Reynolds  number  of 
the  models  is  sufficiently  large. 


SECTION  4 


THE  EXPERIMENTAL  SET-UP  AND  PROCEDURES 

The  flow  induced  by  a  large  circular  area  fire  was  simulated  in  the 
present  experiments  using  two  circular  area  sources  of  helium  mixed  with 
0.5  percent  of  methane  as  a  tracer  gas.  The  combined  relative  density 
of  the  gas  with  respect  to  air  was  p  /p  =0.15.  The  models  were 

O  3 

placed  in  a  5. 5m  x  7. 8m  x  2. 5m  room  in  the  Fluid  Dynamics  and 
Diffusion  Laboratory  (FDDL)  at  Colorado  State  University,  as  shown 
schematically  in  Fig.  3. 

The  pressure  loss  across  the  porous  plate  was  many  times  larger 
than  the  expected  pressure  difference  along  the  plate  due  to  the  induced 
flow,  and  thus  a  fairly  uniform  flux  of  helium  per  unit  area  was 

obtained.  Measurements  with  air  showed  that  the  standard  deviation  of 
the  exit  velocities  over  the  plate  from  the  average  exit  velocity  wq, 
caused  by  the  nonhomogeneity  of  the  plate,  was  of  the  order  of 

8.5  percent. 

The  radius  of  the  plate  in  the  small  circular  model  was 

R  =  26.7  cm  and  the  exit  velocity  of  the  helium  mixture  for 
m 

Q  =  9  SL/s  was  w  =  0.04  m/s,  giving  a  buoyancy  flux  per  unit  area 

2  3  0 

<J  =  0.335  m  /s  .  The  exit  velocity  in  the  large  model  was  even 

smaller. 

Measurements  of  induced  radial  velocities  around  the  fire  were  made 
with  a  calibrated  hot-wire  anemometer.  The  output  of  the  anemometer  was 
fed  directly  to  a  Hewlett-Packard  Series  1000  Real  Time  Executive  Data 
Acquisition  System  of  the  FDDL. 

The  dilution  of  the  helium  mixture  at  different  locations  above  the 
porous  plate  was  determined  by  measuring  the  dilution  of  the  methane 
using  a  gas  chromatograph  with  a  flame  ionization  detector.  The 
relative  density  in  the  model  is  related  to  the  concentration  dilution 
by  (Ap/pa)m=  (ApQ/Pa)  C/Cq.  Samples  of  the  mixture  were  withdrawn  from 
desired  locations  using  1  mm  ID  tubes.  The  sampling  time  was  30  s  and 
thus  the  measured  concentrations  correspond  to  average  values  at  each 
location. 

The  difference  between  the  density  of  the  plume  and  the  density  of 
the  ambient  air  made  it  possible  to  observe  the  shape  of  the  plume  by 
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the  shadowgraph  method.  A  beam  of  light  was  emitted  from  a  projector 
outside  the  room.  The  light  entered  the  room  through  a  window  and  was 
reflected  by  a  large  mirror,  crossed  the  space  above  the  porous  plate 
and  produced  a  shadowgraph  on  a  vertical  screen.  A  wire  grid  placed 
close  to  the  porous  plate  made  it  possible  to  simultaneously  project  a 
grid  on  the  same  screen.  The  angle  of  the  light  beam  was  approximately 
3°  and  thus  the  distortion  due  to  the  distances  between  the  fire,  the 
grid,  and  the  screen  was  small. 

The  refraction  of  light  by  density  gradients  at  the  edges  of  the 
plume  produced  a  dark  shadow  of  the  plume.  Still  photographs  and  a 
video  tape  were  made  of  the  shadowgraph,  from  which  the  visual 
boundaries  of  the  plume,  bs,  were  determined. 

Most  of  the  data  presented  in  the  report  has  been  collected  in  the 
R  =  0.267  m  circular  model.  To  increase  the  induced  radial  velocities 

and  the  Reynolds  number  of  the  model,  large  values  of  buoyancy  flux  has 

2  3 

been  used  in  most  of  the  experiments  in  this  model  (a  =  0.335  m  /s  ). 

Some  measurements  have  been  made  in  a  larger  circular  model  with 

R  =  0.5  m  and  a  buoyancy  flux  of  CT^  =  0.036  m  /s  ).  The  buoyancy  flux 

produced  by  large  urban  fite  has  been  estimated  to  be  of  the  order  of 
2 

(Jp  =  2.1  m  /s.  Thus,  the  second  model  satisfies  the  requirement 
^Ap/pp  =  *  f°r  PrototyPe  fires  with  R  =  1700  m.  Unfortunately,  the 
size  of  the  large  circular  model  was  not  sufficiently  small  relative  to 
the  size  of  the  experimental  chamber,  for  an  accurate  simulation  of  an 
unconfined  burning  fire.  It  is  thus  suspected  that  some  of  the  measure¬ 
ments,  and  in  particular  the  induced  radial  velocities  outside  the 
plume,  could  have  been  affected  by  the  induced  circulation  in  the  room. 

To  get  a  better  picture  of  the  structure  of  the  radial  boundary 
layer,  a  two-dimensional  fire  model  was  built  in  the  test  section  of  the 
Meteorological  Wind  Tunnel  at  the  FDDL.  A  schematic  description  of  the 
experimental  configuration  is  presented  in  Fig.  4. 

The  flow  in  the  outer  region  of  this  two-dimensional  model  is 
expected  to  be  similar  to  the  flow  at  the  outer  edge  of  the  radial 
boundary  layer.  Namely,  Eqs.  ( 23 ) — (28 )  should  apply  to  both  cases. 

To  produce  a  constant  heat  flux,  a  constant  power  per  unit  area  was 


supplied  to  the  heating  coils  below  the  aluminum  plate.  The  large  heat 
flux  into  the  tunnel  made  it  necessary  to  cool  the  air  at  the  two  ends 
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of  the  closed  test  section  of  the  wind  tunnel,  to  avoid  a  rapid 
temperature  rise  in  the  tunnel.  Thus,  the  induced  velocity  toward  the 
center  of  the  plume  in  this  experimental  configuration  is  not  neces¬ 
sarily  equal  to  that  of  an  unbounded  two-dimensional  fire.  The 
structure  of  the  flow  near  the  edge  of  the  hot  plate  should  be  similar, 
however,  to  the  structure  of  the  flow  in  the  radial  boundary  layer  of  a 
large  circular  fire. 


SECTION  5 


THE  EXPERIMENTAL  RESULTS 


5.1  THE  SHAPE  OF  THE  FIRE  PLUME. 

A  photograph  of  the  plume  shadowgraph  is  presented  in  Fig.  5.  The 


average  boundaries  of  the  plume  in  the  shadowgraph,  bs,  are  shown  in 


Fig.  6.  Unless  otherwise  mentioned ,  the  data  shown  in  the  report  has 
been  measured  in  the  small  model  (R  =  0.267  m) !  The  figure  also  shows 
the  radial  distribution  of  the  relative  concentrations  of  the  helium 
mixture  at  several  heights  above  the  source  for  the  small  model.  The 
local  physical  conditions  corresponding  to  the  boundaries  of  the  plume 
in  the  shadowgraph  were  not  established,  but  it  appears  from  the  figure 
that  they  fall  approximately  near  the  edge  of  the  plume,  as  defined  by 
the  concentration  profiles. 

An  attempt  was  made  to  determine  the  plume  width  by  continuously 
traversing,  at  a  low  speed,  a  hot-wire  probe  across  the  plume  and 
recording  the  hot-wire  output  versus  the  radial  distance.  A  typical 
result  of  such  a  recording  is  presented  in  Fig.  7.  The  figure  clearly 
exhibits  strong  turbulence  in  the  convective  plume,  which  has  also  been 
observed  in  the  shadowgraphs.  The  radius  of  the  plume,  ol  ained  by  this 
method,  was  larger  than  the  radius  obtained  by  the  shadowgraph 
technique.  One  possible  reason  for  this  difference  is  the  meandering  of 
the  instantaneous  plume  around  the  average  plume  axis.  These  recordings 


also  suggest  that  while  the  mean  concentrations  and  velocities  across 
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the  convective  plume  are  distributed  as  exp(-r  /2a  ),  the  concentrations 


and  velocities  within  the  instantaneous  plume  are  more  uniformly 
distributed. 

The  flow  visualization  clearly  shows  that  the  plume  of  the  fire 
model  pinched  considerably  due  to  the  induced  radial  velocities,  before 
rising  up.  The  minimum  radius  of  the  visual  plume  in  the  shadowgraph 
was  at  approximately  z  =  17.5  cm,  which  gives 


Z  /R  ~  2/3 
P 


(30) 


The  value  of  the  minimum  radius  of  the  visual  plume  was  approximately 
11.5  cm  or 


b  .  /R  =  0.43 
s  min 


(31) 
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This  observed  contraction  of  the  plume  from  an  area  source  is  consistent 
with  the  visualization  of  the  plume  above  a  model  of  a  heat  island  in 
water  by  Faust  (1981),  although  the  flow  in  that  model  was  partially 
laminar,  with  observations  of  the  fire  plume  in  Project  Flambeau 
(Palmer,  1981),  see  Fig.  8,  and  with  the  theoretical  considerations  of 
Carrier  et  al.  (1984). 


The  concentration  profile  measured  at  z  =  31.7  cm,  z/R  ~  1 . 2  cm 


above  the  porous  plate,  for  a  =  0.335  m  /s  had  a  maximum  concentra¬ 


tion  of  (C/C  )  =  0.375  and  b.  c  =  7.5  cm.  It  follows  that  at  this 

o  max  0.5 


height 


o'  =0^=0- - 2  =  3  13  m/s' 

8max  8p  8C  p  J  J  ' 
a  o  ra 


Since  according  to  Eqs.  (4)  and  (8)  z~" J  =  ll(rtaR^)^^/g'  ,  it 

mdX  max 

follows  that  this  position  corresponds  to  z  =  75.4  cm,  which  is  consis¬ 
tent  with  our  measurement  of  b^  ^  =  7.5  cm,  as  according  to  Eq.  (6) 


bo  5  =  2*  0ne  may  conclude  that  the  virtual  source  of  the  convec¬ 

tive  plume  generated  above  the  fire  model  is  below  the  fire  at 
approximately 


Z  /R  =  -1.64 
v 


The  measurements  in  the  large  model  gave  similar  results,  but  pinching 


was  slightly  larger  and  the  value  of  Zy/R  was  estimated  to  be  in  the 


range  1.4  -  1.6.  It  is  suspected,  however,  that  the  shape  of  the  plume 
in  the  large  model  was  affected  by  the  reduced  circulation  in  the  room. 


5.2  INDUCED  RADIAL  VELOCITIES. 


Induced  radial  velocities  in  the  circular  models  were  measured  at 


r  =  R.  The  magnitude  of  these  velocities  is  very  small,  of  the  order  of 
10  cm/s.  Thus,  the  accuracy  of  these  measurements  is  only  fair, 
±30  percent,  primarily  due  to  background  convective  currents  introduced 
by  nonuniform  temperatures.  It  must  also  be  realized  that  the  radial 
boundary  layer  flow  which  develops  outside  the  fire  region,  is  not 
properly  scaled  in  the  model.  Thus,  these  measurements  can  give  only  a 
rough  estimate  of  the  induced  velocities  in  a  large-area  fire. 
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The  maximum  measured  radial  velocity  outside  the  fire  model  at 
r  =  R,  will  be  referred  to  as  the  induced  radial  velocity  U^.  Figure  9 
shows  the  measured  variation  of  with  the  buoyancy  flux  a  in  the 
small  model,  which  is  closely  described  by 


U. 


R 


(aR) 


1/3 


=  k 


(35) 


with  k  =  0.24.  It  is  noted  that  if  one  assumes  that  the  velocity  field 
above  a  certain  height,  say  Z  ,  is  described  by  Eq.  (3 ) - (8 )  with  a 

3 

virtual  origin  at  Z^,  and  that  the  contribution  of  the  helium  or  burned 
fuel  to  the  mass  balance  is  negligible,  the  average  value  of  the  induced 
velocities  U_  through  the  ring  r  =  R  and  z  =  0  to  Z  can  Be  calculated 
from  the  equation 

U  (Z  /R)/(aR)1/3  =  0.036  [(Za  -  Z  )/R]5/3/(Z  /R)  .  (36) 
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Such  a  calculation  gives  for  T.J R  =  0.75  and  Zy/R  =  -1.5  a  value  of 
0.18,  which  is  consistent  with  our  measurements  as  shown  in  Fig.  10. 
The  calculation  also  indicates  that  models  which  predict  negligible 
pinching,  or  large  negative  values  of  Z^/R,  such  as  the  model  of 
Priestley  and  Ball  (1955)  and  the  numerical  models  of  Small  et  al. 
(1980,  1984),  would  also  predict  higher  induced  radial  velocities. 

Figure  11  shows  the  dimensionless  radial  velocities,  -u  (R,z )/ 

1/3  r 

(ctR)  ,  measured  in  the  two  models.  Clearly  the  values  of  the  dimen¬ 

sionless  radial  velocities  in  the  region  z/R  <  0.3  are  larger  in  the 
large  model.  It  is  suspected  that  the  difference  between  the  two 
models  is  primarily  due  to  induced  circulation  in  the  experimental 
chamber  during  the  experiments  with  the  large  model.  Another  cause  for 
this  difference  is  the  lack  of  similarity  in  the  modeling  of  the 
tuihu.ent  boundary  layer  outside  at  r  >  R. 

In  spite  of  the  differences  between  the  two  models,  it  is  proposed 
that  Eq.  (36)  provides  a  correct  estimate  for  the  order  of  magnitude  of 
the  maximum  induced  radial  velocities  at  r  =  R,  except  that  k  might 
be  larger  by  a  factor  of  two,  namely 

0.24  <  k  <  0.5  .  (37) 


The  figure  also  shows  the  dimensionless  velocities  measured  in  the 
Flambeau  fire  and  its  model  by  Parker  et  al.  (1968),  which  appear  to  be 


consistent  with  the  above  estimate  of  k.  On  the  other  hand,  the  data  in 
the  numerical  model  of  Small  and  Larson  (1984),  which  predicts  only 
small  pinching,  gives  a  value  of  k  =  1.7. 

According  to  Carrier  et  al.  (1985),  the  appropriate  value  of  cr 
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for  the  Hamburg  fire  is  approximately  o  =  2.1  m  /s  .  (o  is  related  to 
e,  the  volumetric  flux  equivalent  of  the  fire  strength  in  that  paper,  by 
O  =  ge/AREA.)  Using  this  value  of  o,  our  estimate  for  full-scale  fires 
with  R  =  2  km  in  neutrally  stratified  atmosphere  gives  radial  wind 
speeds  4  m/s  -8  m/s,  much  lower  than  the  wind  speeds  which  have  been 
associated  with  such  fires. 

5.3  THE  RADIAL  BOUNDARY  LAYER. 

As  seen  in  Fig.  6,  the  growth  of  the  radial  boundary  layer  with 
x  =  R-r  is  rather  fast.  The  angle  Of  in  Fig.  6  is  close  to  34°.  This 
result  appears  to  be  consistent  with  the  shape  of  the  fire  in  Project 
Flambeau,  shown  in  Fig.  8.  As  mentioned  earlier,  the  physical  meaning 
of  the  observed  boundaries  of  the  plume  in  either  the  shadowgraph  method 
or  in  photographs  of  the  full-scale  fire  is  not  very  clear  and  thus  one 
cannot  relate  the  thickness  of  the  radial  boundary  layer  in  Fig.  6  to 
6,  as  defined  in  Eq.  (27). 

One  of  the  questions  which  we  have  tried  to  answer  in  the 
investigation  concerns  the  nature  of  the  radial  boundary  layer,  the 
mixing  within  the  layer,  and  the  availability  of  oxygen  to  the  combus¬ 
tion  layer.  For  this  purpose  one  has  to  estimate  the  relative  signifi¬ 
cance  of  the  turbulence  generated  by  the  shear  stresses  and  the  turbu¬ 
lence  generated  by  the  buoyant  heat  flux.  The  characteristic  scale  of 
the  turbulent  velocities  generated  by  the  buoyancy  flux  is 

tf*(z)  =  (ctz )1/3  ,  (38) 

whereas  the  shear  generated  turbulence  is  of  the  order  of  the  shear 
velocity,  v*,  which  for  a  rough  urban  surface  is  of  the  order  of  U^/10. 
Using  Eq.  (35)  with  k  =  0.4  one  finds  that 

S  25  (|)1/3  .  (39) 

Thus,  except  very  close  to  the  ground  the  buoyancy  generated  turbulence 
will  be  considerably  larger  than  the  shear  generated  turbulence.  In 


other  words,  one  may  positively  conclude  that  this  radial  boundary  layer 
is  a  strong  convective  boundary  layer  (CBL),  similar  to  the  atmospheric 
CBL,  in  which  shear  stresses  and  horizontal  advection  by  U.,  have  an 
insignificant  effect  on  the  mixing  process. 

On  the  basis  of  previous  investigations  of  the  atmospheric  CBL,  the 
heat  and  buoyancy  flux  from  the  combustion  layer  of  a  large  fire  is 
expected  to  be  primarily  in  the  form  of  large  thermals,  rising  at  a 
velocity  slightly  larger  than  w*.  [The  reader  is  referred  to  extensive 
reviews  of  the  CBL  and  diffusion  in  the  CBL  by  Caughey  (1981)  and  Lamb 
(1981)  in  Nieuwstadt  and  von  Dop  (1982).] 

To  compensate  for  the  upward  mass  flux  by  the  thermals,  the 
surrounding  cooler  air  will  descend  at  a  relatively  low  speed.  The 
random  distribution  of  thermals  and  the  descending  air  mass  will  mani¬ 
fest  itself,  at  a  given  point,  by  a  skewed  probability  density  distribu¬ 
tion  (pdd)  of  temperature  (or  density)  fluctuations  and  of  turbulent 
velocity  fluctuations,  which  will  have  a  negative  mode.  This  structure 
creates  a  unique  mass  transfer  pattern,  which  cannot  be  described  by 
gradient-type  numerical  models.  For  example,  plumes  from  elevated 
sources  will  descend  rapidly  toward  the  ground  (as  in  cases  of  fumiga¬ 
tion)  whereas  plumes  from  ground-level  sources  will  be  entrained  into 
thermals  and  rise,  so  that  the  trajectories  of  plumes  from  elevated  and 
from  ground-level  sources  will  appear,  from  the  side,  to  cross  each 
other  (Willis  and  Deardorff,  1976,  1978,  1981;  Poreh  and  Cermak,  1984, 
1985). 

Flow  visualization  and  measurements  in  the  fire  models  have  indeed 
confirmed  the  convective  nature  of  the  flow  in  the  outer  region  of  the 
fire.  Figure  12  shows  a  typical  recording  of  the  temperature  fluctu¬ 
ations  with  time,  made  in  the  outer  region  of  the  two-dimensional  fire 
model.  The  measured  pdd  of  T'/Tj^g  is  shown  in  Fig.  13.  The  shapes  of 
the  fluctuating  signal  and  of  the  dimensionless  pdd  are  similar  to  those 
measured  in  a  simulated  CBL  under  an  elevated  inversion  (Poreh  and 
Cermak,  1985).  They  show  large  positive  spikes  of  T' ,  of  a  relatively 
short  duration,  and  subsidence  of  cooler  air  during  60  percent  of  the 
time . 

Visualization  of  the  flow  in  the  outer  region  of  the  two- 
dimensional  model  was  achieved  by  spreading  titanium  tetrachloride  (TTC) 


on  the  heated  aluminum  plate.  The  existence  and  motion  of  thermals  were 
clearly  observed  and  were  recorded  on  video  tape.  The  spherical  shape 
of  the  rising  thermals  can  also  be  observed  in  the  still  photographs 
presented  in  Figs.  14  and  15.  (The  flow  in  all  the  figures  is  frcm 
right  to  left.)  Figure  15b  shows  the  central  region  of  the  fire  plume. 
Since  the  smoke  was  introduced  on  the  right-hand  side  of  the  model,  the 
flow  on  the  left  of  the  centerline  is  not  observed. 

Figure  15b  shows  smoke  emitted  from  ground  level  in  the  outer 
region  together  with  a  smoke  filament  emitted  from  an  elevated  source 
above  the  convective  layer.  If  the  developing  boundary  layer  below  the 
smoke  filament  were  shear  dominated,  the  smoke  filament  would  have  been 
rapidly  broken,  upon  entering  the  developing  layer,  but  its  center  of 
gravity  would  have  not  significantly  changed.  The  photographs  presented 
in  Fig.  17,  taken  after  the  TTC  spread  on  the  floor  bed  evaporated,  show 
what  happened,  most  of  the  time  to  smoke  filaments  released  above  the 
radial  boundary  layer.  When  they  intersected  the  developing  boundary 
layer,  they  were  carried  down  toward  the  hot  floor,  and  only  later 
entrained  into  the  hot  thermals  and  carried  up.  In  many  cases  such 
filaments  reached  the  center  core  of  the  fire  model  before  rising, 
although  many  times  smoke  filaments  were  entrained  into  a  rising  thermal 
before  reaching  the  center  core.  The  same  pattern  of  fresh  air 
entrainment  was  observed  in  the  circular  fire  models. 


To  examine  the  effect  of  a  heterogeneous,  cellular  distribution  of 
fuel  on  the  plume  shape,  the  porous  plate  in  the  small  circular  model 
was  covered  by  large  plates  with  large  perforations:  (1)  12  percent 
porosity,  1  cm  diameter  holes  and  (2)  30  percent  porosity,  1.5  cm 
diameter  holes.  One  could  not  recognize  changes  in  the  shape  of  the 
plume  due  to  the  heterogeneous  supply  of  helium.  We  shall  later  see 
that  similar  results  were  obtained  in  small  fires. 

Could  swirling  of  the  flow  in  the  radial  boundary  layer,  induced 
by  a  preexisting  circulation  in  the  region,  change  the  convective  struc¬ 
ture  of  the  outer  radial  layer?  The  shear  induced  by  a  circulatory 
motion  will  start  to  dominate  the  mixing  process  in  the  radial  layer 
only  if  the  tangential  velocity  will  be  at  least  one  o  der  of 

magnitude  larger  than  U^.  In  such  cases,  centripetal  acceleration 
acting  on  the  relatively  lighter  thermals,  which  is  proportional  to 
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ApV^/ (pr) ,  will  increase  and  become  significant  with  respect  to  the 
vertical  acceleration  g' .  Obviously,  such  an  effect  will  be  more 
significant  near  the  core  of  the  fire  and  in  the  convective  plume  where 
r  is  small  than  in  the  radial  boundary  layer. 


5.4  DILUTION  OF  THE  BUOYANT  GASES. 

Concentration  measurements,  made  along  the  axis  of  symmetry  of  the 
plumes  in  the  two  circular  models,  were  used  to  calculate  the  distribu¬ 
tion  of  the  dimensionless  buoyancy  factor  along  the  vertical  axis  of  the 
model  (see  Fig.  18).  The  figure  also  shows  the  values  of  Bq  which 
correspond  to  the  source  gas  in  the  two  models. 

The  measurements  indicate  that  very  rapid  mixing  of  the  buoyant 
gases  occurs  near  the  ground.  The  relative  thickness  of  pseudo¬ 
combustion  layer  in  the  model  is  very  small,  of  the  order  of  0.03  in  the 
small  model  and  smaller  than  0.005  in  the  large  model.  The  measurements 
are  also  consistent  with  temperature  measurements  in  the  two-dimensional 
model.  Although  the  temperature  difference  between  the  hot  floor  and 
the  ambient  air,  in  this  model,  was  of  the  order  of  130°C,  the  tempera¬ 
ture  increase  of  the  air  at  a  height  of  0.5  cm  above  the  hot  floor  at 
x  =  50  cm  from  the  edge  of  the  hot  plate  was  only  18°C. 

The  figure  also  shows  the  dimensionless  buoyancy  factor  in  small 
burning  alcohol  fires  (Yokoi,  1961).  The  relevance  of  Yokoi's  study  to 
our  problem  was  not  recognized  earlier  due  to  the  title  of  his  paper. 
His  experiments  may  be  considered,  however,  to  be  models  of  LAF,  except 
that  the  relative  thickness  of  the  combustion  layer  in  his  models  was 
probably  larger  than  in  LAF  and  some  of  the  energy  Qg  was  lost  by 
radiation. 

In  spite  of  these  shortcomings,  the  distribution  of  the 
dimensionless  buoyancy  factor  B  (denoted  in  Yokoi's  paper  by  (h)  )  is 
very  similar  to  that  measured  in  the  present  study.  Both  studies  sug¬ 
gest  that  B  is  approximately  constant  in  the  core  of  the  fire  above  the 
combustion  layer.  The  slightly  smaller  values  of  B  in  Yokoi's  models 
is  probably  due  to  radiation  losses.  When  analyzing  the  far  field  data, 
Yokoi  assumed  that  the  virtual  origin  of  the  convective  plume  is  at 
Z  =  0.  As  shown  in  the  figure,  his  data  match  much  better  Eqs.  (2)  and 
(6)  with  Zv/R  S  1. 


Yokoi  has  also  studied  a  circular  area  source  made  out  of  many 
little  wicks  of  alcohol  lamps,  which  he  termed  "discontinuous  heat 
source".  As  the  data  in  this  figure  indicate,  no  difference  was 
observed  between  the  continuous  and  discontinuous  sources.  (We  did  not 
show  in  the  figure  his  measurements  in  a  R  =  3.3  cm  fire  and  a  R  =  6  cm 
fire,  which  we  observed  to  be  partially  laminar.)  Yokoi  also  modeled 
rectangular  fires,  which  gave  slightly  larger  initial  mixing. 

The  approximate  equality  of  the  dimensionless  buoyancy  factors  in 
these  studies  supports  our  basic  assumption  that  the  induced  velocities 
and  the  buoyancy  factor  gAp/p  above  the  combustion  layer  are  functions 
of  the  buoyancy  flux  a  and  the  size  of  the  large  fire  R,  and  inde¬ 
pendent  of  the  initial  density  of  the  buoyant  gases.  Note  that  the 
initial  dimensionless  buoyancy  in  the  two  helium  models  differed  by  a 
factor  of  five. 

Since  g'  is  approximately  proportional  to  AT/T  it  follows  that 

-1/3 

AT/T  in  large-area  fire  plumes  decreases  as  R  '  .  Namely,  the  aver¬ 
age  temperature  rise,  as  well  as  the  depletion  of  oxygen  decrease  with 
the  size  of  the  fire.  It  should  be  stressed  that  this  conclusion 
applies  only  to  average  values.  The  local  temperature  rise  or  oxygen 
depletion  in  an  individual  thermal,  or  in  an  individual  fire  plume, 
would  be  governed  by  the  individual  fire  size  and  rate  of  burning. 

The  order  of  magnitude  of  the  average  excess  temperatures  in  the 
central  core  of  LAF  has  been  calculated  and  plotted  in  Fig.  19  assuming 
that  the  values  of  B  above  the  combustion  layer  (z/R  =  0.005)  and 

above  the  pinch  point  are  3.2  and  3.3  (see  Fig.  18),  and  that  the  value 
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of  the  buoyancy  flux  is  a  =  2.1  m  /s  .  The  estimate  clearly  shows  that 
the  average  temperature  rise  in  large-area  fire  plumes  is  not  going  to 
be  large. 

Again,  the  results  are  initially  surprising,  as  one  intuitively 
expects  the  excess  temperatures  to  increase  with  the  size  of  the  fire. 
This  intuitive  correlation  with  the  size  of  the  fire  is,  however,  mis¬ 
leading,  as  the  temperature  rise  is  determined  by  the  ratio  of  the  flux 

of  burned  gases  to  the  flux  of  fresh  air.  The  flux  of  burned  gases  is 

2 

proportional  to  aR  whereas  the  flux  of  fresh  air  is  proportional  to 
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U_R  .  Since  UD  is  proportional  to  (aR)  ,  the  ratio  of  these  fluxes 

K  K  -1/3 

is  proportional  to  R  . 


5.5  ESTIMATE  OF  PLUME  RISE  FROM  AREA  SOURCES . 


The  fact  that  the  large-area  fires  have  convective  plumes  that 
behave  at  large  heights  as  weak  convective  plumes  from  a  virtual  origin 
at  z  =  -Zv,  suggest  that  the  rise  of  such  a  plume  in  a  stably  stratified 
atmosphere  would  be  much  smaller  than  that  of  a  plume  with  an  equal 
total  buoyancy  flux  rising  from  a  point  source.  An  approximate  estimate 
of  the  maximum  rise  of  a  plume  from  an  area  source  can  be  made  by  assum¬ 
ing  that  at  z  =  0  the  width  of  the  plume  is  equal  to  that  of  a  plume 

which  originates  at  the  virtual  origin  Z  and  which  at  z  =  Z  (Z  =  0) 
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has  a  momentum  flux  tojR  .  The  value  of  z,  see  Eq.  (8),  corresponding 

to  ground  level  for  such  a  plume,  can  be  calculated  from 


'1 


f1/4  (z:) 


a 1,2  Z 

-3/8 

G  v 

(noR2)1/4 


(40) 


where  f(z^)  is  the  dimensionless  flux  defined  in  Eq.  (7). 

Such  a  model  suggests  a  drastic,  and  initially  surprising, 

reduction  of  the  maximum  plume  rise.  Figure  20  shows  the  maximum  plume 

-4  2 

rise  as  a  function  of  Z  ,  for  G  =  1.2  x  10  1/s  (approximately 

V  2 

dO/dZ  =  0.0035°C/m),  for  a  100  km  fire.  (Baseline  data,  as  given  in 
the  National  Research  Council  Report  (1985).  A  maximum  plume  rise  of 
10  km  is  proposed  in  the  report  for  such  a  fire.) 

The  initial  surprise  diminishes  when  one  recalls  that  the 
theoretical  solution  of  Morton  et  al.  (1956)  indicates  that  most  of  the 
buoyancy  flux  of  a  plume  rising  in  a  stratified  atmosphere  occurs  at 
Zj  >  1,  after  the  plume  has  acquired  a  large  width.  The  initial  finite 
size  of  the  plume  from  a  large-area  fire  is  therefore  the  primary  cause 
for  the  reduction  of  the  maximum  plume  rise. 

An  interesting  corollary  is  that  models  which  predict  large  induced 
radial  velocities  (a  large  value  of  -Zv),  should  also  predict  a  small 
plume  rise,  and  vise-versa! 

5.6  VIDEOTAPE  RECORDINGS  OF  THE  FLOW  VISUALIZATION. 

Flow  visualization  experiments  of  the  flow  in  the  convective  layer 
of  the  two-dimensional  models  and  of  the  shape  of  the  plume  for  dif¬ 
ferent  configurations  of  area  sources  were  recorded  on  video  tape  and 
are  discussed  in  a  separate  report  (Stout,  1986). 
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SECTION  6 


CONCLUSIONS 

Physical  simulations  of  plumes  from  circular,  large-area  fires, 
using  a  0.53  m  diameter  and  a  1.00  m  diameter  area  sources  of  helium, 
and  a  simulation  of  a  two-dimensional  fire  using  a  2m  x  2m  hot  plate, 
have  yielded  important  information  about  the  characteristics  of  such 
plumes  and  the  nature  of  the  flow  above  the  combustion  layer. 

Regions  of  distinct  flow  characteristics  were  identified:  a 
relatively  shallow  combustion  layer  near  the  ground,  a  radial  convective 
boundary  layer,  a  central  core ,  and  a  convective  plume  (see  Fig.  1).  It 
was  shown  that  large-area  fires  induce  radially  inward  winds  in  the 
surface  layer  around  the  fire,  which  are  proportional  to  the  cube  root 
of  the  radius  of  the  fire  and  of  the  buoyancy  flux.  The  magnitude  of 
the  induced  velocities  near  the  ground  are  larger  than  those  induced  by 
point  sources  with  the  same  buoyancy  flux,  however  they  are  small 
compared  to  previous  estimates. 

The  structure  of  the  radial  convective  boundary  layer,  which 
develops  in  the  outer  region  of  large-area  fires  above  the  shallow 
combustion  layer,  has  been  shown  to  be  similar  to  the  convective 
boundary  layer  in  the  atmosphere.  Heat  or  buoyancy  transfer  from  the 
combustion  layer  is  primarily  in  the  form  of  thermals  that  induce  subsi¬ 
dence  of  fresh  air  with  oxygen  toward  the  ground  and  toward  the  central 
core  of  the  fire.  These  characteristics  are  inherent  properties  of  the 
convective  flow  and  should  exist  even  for  an  inhomogeneous,  cellular 
combustion  layer. 

The  radial  flow  causes  a  considerable  contraction  of  the  fire  plume 
before  it  rises.  In  addition,  the  intense  mixing  of  the  gases,  prior  to 
the  formation  of  the  convective  plume,  reduces  drastically  the  excess 
temperatures  at  the  base  of  the  convective  plume,  whose  characteristics 
at  higher  elevations  appear  to  be  similar  to  those  of  weakly  buoyant 
plumes  rising  from  a  virtual  point  source  at  approximately 
z/R  =  0(-l  .5) . 

The  excess  temperatures  in  the  central  core  of  LAF,  above  the 

combustion  layer,  are  found  to  be  approximately  constant  and  propor- 
-1/3 

tional  to  R  .  This  finding  implies  that  the  average  temperature 


rise  and  oxygen  depletion  in  plumes  of  LAF  decrease  with  the  size  of  the 
fire. 


Prior  mixing  in  the  radial  convective  boundary  layer  is  also 
expected  to  decrease  the  rise  of  the  convective  plume  in  upper  stably 
stratified  layers,  below  the  rise  predicted  for  convective  plumes  with 
the  same  total  buoyancy  flux  from  a  point  source  at  ground  level. 

The  results  also  suggest  that  the  burning  process  of  small  fires, 
which  burn  simultaneously  over,  a  large  urban  area,  is  similar  to  that  of 
a  small  single  fire  in  an  atmosphere  with  relatively  small  ambient 
horizontal  velocities.  Thus,  one  may  study  the  smoke  production  pro¬ 
cesses  in  LAF  in  a  relatively  small  fire.  This  suggestion  is  based  on 
the  conclusions  that  the  induced  horizontal  velocities  in  the  large-area 
fires  are  small  compared  to  the  convective  vertical  velocities  of  the 
plumes  from  individual  fires  and  that  the  oxygen  level  in  the  ambient 
air  is  not  drastically  reduced. 

Finally,  the  results  obtained  indicate  that  significant  aspects  of 
large-area  fires  can  be  simulated  conveniently  in  small-scale  physical 
models.  It  is  thus  proposed  to  use  the  same  approach  to  study  the 
effect  of  ambient  cross  winds  and  of  atmospheric  circulation  and  strati¬ 
fication  on  the  flow  field  induced  by  large-area  fires.  It  is  also 
proposed  to  study  the  temporal  and  spatial  variations  of  the  temperature 
field  in  large-area  fires  and  the  detailed  effect  of  the  fuel  bed 
nonhomogeneity  on  the  fire  plume  using  the  same  method. 
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Figure  1.  Schematic  description  of  a  large  fire. 


J.F.R.C..  N.8.S.  AND  U.S. 
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hot-wire  probe  traversing  across  the  tire  pi  nine  (■/. 


Figure  8.  A  15-acre  fire  in  the  Flambeau  Project.  Trees  in 
foreground  are  12  ft  high. 
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Estimate  of  Entraiment  Velocity 
into  the  Convective  Plume 


o  -  Data  (  a  *  0.335  m2/s3 ) 


Estimate  for  Average 
Speed  at 

0  <  z/R  <  3/4. 


Ur  (  R  )  /  (  cr  R  )  ; 


Figure  10.  Radial  velocities  at  the  edge  of  the  small  circular  model 
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Figure  16.  Photographs  of  the  CBL  in  the  two-dimensional  fire  model 
together  with  a  smoke  filament  released  above  the  layer 


Figure  17a. 


Photograph  of  a  smoke  filament  released  above 
in  the  two-dimensional  fire  model. 
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Figure  18.  Vertical  distribution  of  the  dimensionless  buoyancy 
factor  B  along  the  axis  of  different  fire  models. 
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APPENDIX 


LIST  OF  SYMBOLS 


Definition 


Typical  radial  size  of  the  convective  plume. 

Boundaries  of  the  plume  as  determined  by  the  shadowgraph 
method. 

Defined  in  Eq.  (6). 

Dimensionless  buoyancy  factor,  defined  in  Eq.  (23). 
Concentration  of  the  tracer  gas. 

Concentration  of  the  tracer  gas  at  the  source. 

Relative  buoyancy  flux;  see  Eq.  (7). 

=  V"- 

00 

Total  buoyancy  flux  at  =  0,  Fq  =  271  f  wg'  rdr. 

o 

Total  buoyancy  flux  at  Zy 
=  gAp/p  ,  buoyancy  factor. 

(g/0)d0/dz. 

Constant  in  Eq.  (36). 

Monin-Obukhov  length. 

In  Fig.  2.  Flame  height/fire  linear  dimension. 

2 

Brunt-Vaisalla  frequency,  N  =  G. 

Radial  distance. 

Temperature  (°K). 

Any  velocity  component. 

Ambient  crosswind  speed. 

Maximum  inward  radial  speed  at  the  edge  of  the  fire. 

Shear  velocity. 

Radial  component  of  the  flow  field  (toward  the  fire). 


Symbol  Definition 


Vertical  velocity  in  the  convective  plume  region. 


Representative  vertical  velocity  in  the  convective  plume. 


Characteristics  velocity  scale  in  convective  boundary  layer. 


Upward  mean  velocity  of  the  helium  at  the  source. 


Distance  from  the  outer  edge  of  the  fire. 


Height  above  the  ground. 


Height  above  virtual  origin  of  a  convective  plume, 

z  =  2  -  Z  . 

v 


Position  of  virtual  source. 


Dimensionless  height,  defined  in  Eq.  (8). 


Plume  rise. 


Height  above  the  ground  of  the  pitch  point,  where  the  radius 
of  the  fire  plume  is  minimum  (see  Fig.  1). 


=  Gp/5. 


=  2-l/2  a. 


Defined  in  Eq.  (3) 


Maximum  upward  velocity  in  the  convective  plume. 


Defined  in  Eq.  (26). 


=  Pa  -  P- 


Potential  temperature. 


Density  in  the  model  or  potential  density  in  the  atmosphere. 


Buoyancy  flux  per  unit  area. 


Angular  velocity. 


Subscrii 


Denotes  value  in  a  model. 


Denotes  value  in  prototype. 


Denotes  ambient  conditions. 


Denotes  initial  conditions  at  the  source. 
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